3.254 \(\int \frac{x^4}{(d+e x^2) (a+c x^4)^2} \, dx\)

Optimal. Leaf size=685 \[ \frac{\left (\sqrt{a} e+3 \sqrt{c} d\right ) \log \left (-\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{a}+\sqrt{c} x^2\right )}{16 \sqrt{2} a^{3/4} c^{3/4} \left (a e^2+c d^2\right )}-\frac{\left (\sqrt{a} e+3 \sqrt{c} d\right ) \log \left (\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{a}+\sqrt{c} x^2\right )}{16 \sqrt{2} a^{3/4} c^{3/4} \left (a e^2+c d^2\right )}+\frac{\left (3 \sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (1-\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt{2} a^{3/4} c^{3/4} \left (a e^2+c d^2\right )}-\frac{\left (3 \sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}+1\right )}{8 \sqrt{2} a^{3/4} c^{3/4} \left (a e^2+c d^2\right )}-\frac{\sqrt [4]{c} d^2 \left (\sqrt{a} e+\sqrt{c} d\right ) \log \left (-\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{a}+\sqrt{c} x^2\right )}{4 \sqrt{2} a^{3/4} \left (a e^2+c d^2\right )^2}+\frac{\sqrt [4]{c} d^2 \left (\sqrt{a} e+\sqrt{c} d\right ) \log \left (\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{a}+\sqrt{c} x^2\right )}{4 \sqrt{2} a^{3/4} \left (a e^2+c d^2\right )^2}-\frac{\sqrt [4]{c} d^2 \left (\sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (1-\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt{2} a^{3/4} \left (a e^2+c d^2\right )^2}+\frac{\sqrt [4]{c} d^2 \left (\sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}+1\right )}{2 \sqrt{2} a^{3/4} \left (a e^2+c d^2\right )^2}-\frac{x \left (d-e x^2\right )}{4 \left (a+c x^4\right ) \left (a e^2+c d^2\right )}+\frac{d^{3/2} e^{3/2} \tan ^{-1}\left (\frac{\sqrt{e} x}{\sqrt{d}}\right )}{\left (a e^2+c d^2\right )^2} \]

[Out]

-(x*(d - e*x^2))/(4*(c*d^2 + a*e^2)*(a + c*x^4)) + (d^(3/2)*e^(3/2)*ArcTan[(Sqrt[e]*x)/Sqrt[d]])/(c*d^2 + a*e^
2)^2 - (c^(1/4)*d^2*(Sqrt[c]*d - Sqrt[a]*e)*ArcTan[1 - (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/(2*Sqrt[2]*a^(3/4)*(c*d^2
 + a*e^2)^2) + ((3*Sqrt[c]*d - Sqrt[a]*e)*ArcTan[1 - (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/(8*Sqrt[2]*a^(3/4)*c^(3/4)*
(c*d^2 + a*e^2)) + (c^(1/4)*d^2*(Sqrt[c]*d - Sqrt[a]*e)*ArcTan[1 + (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/(2*Sqrt[2]*a^
(3/4)*(c*d^2 + a*e^2)^2) - ((3*Sqrt[c]*d - Sqrt[a]*e)*ArcTan[1 + (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/(8*Sqrt[2]*a^(3
/4)*c^(3/4)*(c*d^2 + a*e^2)) - (c^(1/4)*d^2*(Sqrt[c]*d + Sqrt[a]*e)*Log[Sqrt[a] - Sqrt[2]*a^(1/4)*c^(1/4)*x +
Sqrt[c]*x^2])/(4*Sqrt[2]*a^(3/4)*(c*d^2 + a*e^2)^2) + ((3*Sqrt[c]*d + Sqrt[a]*e)*Log[Sqrt[a] - Sqrt[2]*a^(1/4)
*c^(1/4)*x + Sqrt[c]*x^2])/(16*Sqrt[2]*a^(3/4)*c^(3/4)*(c*d^2 + a*e^2)) + (c^(1/4)*d^2*(Sqrt[c]*d + Sqrt[a]*e)
*Log[Sqrt[a] + Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2])/(4*Sqrt[2]*a^(3/4)*(c*d^2 + a*e^2)^2) - ((3*Sqrt[c]*d
 + Sqrt[a]*e)*Log[Sqrt[a] + Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2])/(16*Sqrt[2]*a^(3/4)*c^(3/4)*(c*d^2 + a*e
^2))

________________________________________________________________________________________

Rubi [A]  time = 0.609404, antiderivative size = 685, normalized size of antiderivative = 1., number of steps used = 23, number of rules used = 10, integrand size = 22, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.454, Rules used = {1314, 1179, 1168, 1162, 617, 204, 1165, 628, 1171, 205} \[ \frac{\left (\sqrt{a} e+3 \sqrt{c} d\right ) \log \left (-\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{a}+\sqrt{c} x^2\right )}{16 \sqrt{2} a^{3/4} c^{3/4} \left (a e^2+c d^2\right )}-\frac{\left (\sqrt{a} e+3 \sqrt{c} d\right ) \log \left (\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{a}+\sqrt{c} x^2\right )}{16 \sqrt{2} a^{3/4} c^{3/4} \left (a e^2+c d^2\right )}+\frac{\left (3 \sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (1-\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt{2} a^{3/4} c^{3/4} \left (a e^2+c d^2\right )}-\frac{\left (3 \sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}+1\right )}{8 \sqrt{2} a^{3/4} c^{3/4} \left (a e^2+c d^2\right )}-\frac{\sqrt [4]{c} d^2 \left (\sqrt{a} e+\sqrt{c} d\right ) \log \left (-\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{a}+\sqrt{c} x^2\right )}{4 \sqrt{2} a^{3/4} \left (a e^2+c d^2\right )^2}+\frac{\sqrt [4]{c} d^2 \left (\sqrt{a} e+\sqrt{c} d\right ) \log \left (\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{a}+\sqrt{c} x^2\right )}{4 \sqrt{2} a^{3/4} \left (a e^2+c d^2\right )^2}-\frac{\sqrt [4]{c} d^2 \left (\sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (1-\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt{2} a^{3/4} \left (a e^2+c d^2\right )^2}+\frac{\sqrt [4]{c} d^2 \left (\sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}+1\right )}{2 \sqrt{2} a^{3/4} \left (a e^2+c d^2\right )^2}-\frac{x \left (d-e x^2\right )}{4 \left (a+c x^4\right ) \left (a e^2+c d^2\right )}+\frac{d^{3/2} e^{3/2} \tan ^{-1}\left (\frac{\sqrt{e} x}{\sqrt{d}}\right )}{\left (a e^2+c d^2\right )^2} \]

Antiderivative was successfully verified.

[In]

Int[x^4/((d + e*x^2)*(a + c*x^4)^2),x]

[Out]

-(x*(d - e*x^2))/(4*(c*d^2 + a*e^2)*(a + c*x^4)) + (d^(3/2)*e^(3/2)*ArcTan[(Sqrt[e]*x)/Sqrt[d]])/(c*d^2 + a*e^
2)^2 - (c^(1/4)*d^2*(Sqrt[c]*d - Sqrt[a]*e)*ArcTan[1 - (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/(2*Sqrt[2]*a^(3/4)*(c*d^2
 + a*e^2)^2) + ((3*Sqrt[c]*d - Sqrt[a]*e)*ArcTan[1 - (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/(8*Sqrt[2]*a^(3/4)*c^(3/4)*
(c*d^2 + a*e^2)) + (c^(1/4)*d^2*(Sqrt[c]*d - Sqrt[a]*e)*ArcTan[1 + (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/(2*Sqrt[2]*a^
(3/4)*(c*d^2 + a*e^2)^2) - ((3*Sqrt[c]*d - Sqrt[a]*e)*ArcTan[1 + (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/(8*Sqrt[2]*a^(3
/4)*c^(3/4)*(c*d^2 + a*e^2)) - (c^(1/4)*d^2*(Sqrt[c]*d + Sqrt[a]*e)*Log[Sqrt[a] - Sqrt[2]*a^(1/4)*c^(1/4)*x +
Sqrt[c]*x^2])/(4*Sqrt[2]*a^(3/4)*(c*d^2 + a*e^2)^2) + ((3*Sqrt[c]*d + Sqrt[a]*e)*Log[Sqrt[a] - Sqrt[2]*a^(1/4)
*c^(1/4)*x + Sqrt[c]*x^2])/(16*Sqrt[2]*a^(3/4)*c^(3/4)*(c*d^2 + a*e^2)) + (c^(1/4)*d^2*(Sqrt[c]*d + Sqrt[a]*e)
*Log[Sqrt[a] + Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2])/(4*Sqrt[2]*a^(3/4)*(c*d^2 + a*e^2)^2) - ((3*Sqrt[c]*d
 + Sqrt[a]*e)*Log[Sqrt[a] + Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2])/(16*Sqrt[2]*a^(3/4)*c^(3/4)*(c*d^2 + a*e
^2))

Rule 1314

Int[(((f_.)*(x_))^(m_.)*((a_) + (c_.)*(x_)^4)^(p_))/((d_.) + (e_.)*(x_)^2), x_Symbol] :> -Dist[(a*f^4)/(c*d^2
+ a*e^2), Int[(f*x)^(m - 4)*(d - e*x^2)*(a + c*x^4)^p, x], x] + Dist[(d^2*f^4)/(c*d^2 + a*e^2), Int[((f*x)^(m
- 4)*(a + c*x^4)^(p + 1))/(d + e*x^2), x], x] /; FreeQ[{a, c, d, e, f}, x] && LtQ[p, -1] && GtQ[m, 2]

Rule 1179

Int[((d_) + (e_.)*(x_)^2)*((a_) + (c_.)*(x_)^4)^(p_), x_Symbol] :> -Simp[(x*(d + e*x^2)*(a + c*x^4)^(p + 1))/(
4*a*(p + 1)), x] + Dist[1/(4*a*(p + 1)), Int[Simp[d*(4*p + 5) + e*(4*p + 7)*x^2, x]*(a + c*x^4)^(p + 1), x], x
] /; FreeQ[{a, c, d, e}, x] && NeQ[c*d^2 + a*e^2, 0] && LtQ[p, -1] && IntegerQ[2*p]

Rule 1168

Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[a*c, 2]}, Dist[(d*q + a*e)/(2*a*c),
 Int[(q + c*x^2)/(a + c*x^4), x], x] + Dist[(d*q - a*e)/(2*a*c), Int[(q - c*x^2)/(a + c*x^4), x], x]] /; FreeQ
[{a, c, d, e}, x] && NeQ[c*d^2 + a*e^2, 0] && NeQ[c*d^2 - a*e^2, 0] && NegQ[-(a*c)]

Rule 1162

Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[(2*d)/e, 2]}, Dist[e/(2*c), Int[1/S
imp[d/e + q*x + x^2, x], x], x] + Dist[e/(2*c), Int[1/Simp[d/e - q*x + x^2, x], x], x]] /; FreeQ[{a, c, d, e},
 x] && EqQ[c*d^2 - a*e^2, 0] && PosQ[d*e]

Rule 617

Int[((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> With[{q = 1 - 4*Simplify[(a*c)/b^2]}, Dist[-2/b, Sub
st[Int[1/(q - x^2), x], x, 1 + (2*c*x)/b], x] /; RationalQ[q] && (EqQ[q^2, 1] ||  !RationalQ[b^2 - 4*a*c])] /;
 FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]

Rule 204

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> -Simp[ArcTan[(Rt[-b, 2]*x)/Rt[-a, 2]]/(Rt[-a, 2]*Rt[-b, 2]), x] /
; FreeQ[{a, b}, x] && PosQ[a/b] && (LtQ[a, 0] || LtQ[b, 0])

Rule 1165

Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[(-2*d)/e, 2]}, Dist[e/(2*c*q), Int[
(q - 2*x)/Simp[d/e + q*x - x^2, x], x], x] + Dist[e/(2*c*q), Int[(q + 2*x)/Simp[d/e - q*x - x^2, x], x], x]] /
; FreeQ[{a, c, d, e}, x] && EqQ[c*d^2 - a*e^2, 0] && NegQ[d*e]

Rule 628

Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Simp[(d*Log[RemoveContent[a + b*x +
c*x^2, x]])/b, x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[2*c*d - b*e, 0]

Rule 1171

Int[((d_) + (e_.)*(x_)^2)^(q_)/((a_) + (c_.)*(x_)^4), x_Symbol] :> Int[ExpandIntegrand[(d + e*x^2)^q/(a + c*x^
4), x], x] /; FreeQ[{a, c, d, e}, x] && NeQ[c*d^2 + a*e^2, 0] && IntegerQ[q]

Rule 205

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]*ArcTan[x/Rt[a/b, 2]])/a, x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rubi steps

\begin{align*} \int \frac{x^4}{\left (d+e x^2\right ) \left (a+c x^4\right )^2} \, dx &=-\frac{a \int \frac{d-e x^2}{\left (a+c x^4\right )^2} \, dx}{c d^2+a e^2}+\frac{d^2 \int \frac{1}{\left (d+e x^2\right ) \left (a+c x^4\right )} \, dx}{c d^2+a e^2}\\ &=-\frac{x \left (d-e x^2\right )}{4 \left (c d^2+a e^2\right ) \left (a+c x^4\right )}+\frac{\int \frac{-3 d+e x^2}{a+c x^4} \, dx}{4 \left (c d^2+a e^2\right )}+\frac{d^2 \int \left (\frac{e^2}{\left (c d^2+a e^2\right ) \left (d+e x^2\right )}+\frac{c \left (d-e x^2\right )}{\left (c d^2+a e^2\right ) \left (a+c x^4\right )}\right ) \, dx}{c d^2+a e^2}\\ &=-\frac{x \left (d-e x^2\right )}{4 \left (c d^2+a e^2\right ) \left (a+c x^4\right )}+\frac{\left (c d^2\right ) \int \frac{d-e x^2}{a+c x^4} \, dx}{\left (c d^2+a e^2\right )^2}+\frac{\left (d^2 e^2\right ) \int \frac{1}{d+e x^2} \, dx}{\left (c d^2+a e^2\right )^2}-\frac{\left (\frac{3 \sqrt{c} d}{\sqrt{a}}-e\right ) \int \frac{\sqrt{a} \sqrt{c}+c x^2}{a+c x^4} \, dx}{8 c \left (c d^2+a e^2\right )}-\frac{\left (\frac{3 \sqrt{c} d}{\sqrt{a}}+e\right ) \int \frac{\sqrt{a} \sqrt{c}-c x^2}{a+c x^4} \, dx}{8 c \left (c d^2+a e^2\right )}\\ &=-\frac{x \left (d-e x^2\right )}{4 \left (c d^2+a e^2\right ) \left (a+c x^4\right )}+\frac{d^{3/2} e^{3/2} \tan ^{-1}\left (\frac{\sqrt{e} x}{\sqrt{d}}\right )}{\left (c d^2+a e^2\right )^2}+\frac{\left (d^2 \left (\frac{\sqrt{c} d}{\sqrt{a}}-e\right )\right ) \int \frac{\sqrt{a} \sqrt{c}+c x^2}{a+c x^4} \, dx}{2 \left (c d^2+a e^2\right )^2}+\frac{\left (d^2 \left (\frac{\sqrt{c} d}{\sqrt{a}}+e\right )\right ) \int \frac{\sqrt{a} \sqrt{c}-c x^2}{a+c x^4} \, dx}{2 \left (c d^2+a e^2\right )^2}-\frac{\left (\frac{3 \sqrt{c} d}{\sqrt{a}}-e\right ) \int \frac{1}{\frac{\sqrt{a}}{\sqrt{c}}-\frac{\sqrt{2} \sqrt [4]{a} x}{\sqrt [4]{c}}+x^2} \, dx}{16 c \left (c d^2+a e^2\right )}-\frac{\left (\frac{3 \sqrt{c} d}{\sqrt{a}}-e\right ) \int \frac{1}{\frac{\sqrt{a}}{\sqrt{c}}+\frac{\sqrt{2} \sqrt [4]{a} x}{\sqrt [4]{c}}+x^2} \, dx}{16 c \left (c d^2+a e^2\right )}+\frac{\left (3 \sqrt{c} d+\sqrt{a} e\right ) \int \frac{\frac{\sqrt{2} \sqrt [4]{a}}{\sqrt [4]{c}}+2 x}{-\frac{\sqrt{a}}{\sqrt{c}}-\frac{\sqrt{2} \sqrt [4]{a} x}{\sqrt [4]{c}}-x^2} \, dx}{16 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}+\frac{\left (3 \sqrt{c} d+\sqrt{a} e\right ) \int \frac{\frac{\sqrt{2} \sqrt [4]{a}}{\sqrt [4]{c}}-2 x}{-\frac{\sqrt{a}}{\sqrt{c}}+\frac{\sqrt{2} \sqrt [4]{a} x}{\sqrt [4]{c}}-x^2} \, dx}{16 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}\\ &=-\frac{x \left (d-e x^2\right )}{4 \left (c d^2+a e^2\right ) \left (a+c x^4\right )}+\frac{d^{3/2} e^{3/2} \tan ^{-1}\left (\frac{\sqrt{e} x}{\sqrt{d}}\right )}{\left (c d^2+a e^2\right )^2}+\frac{\left (3 \sqrt{c} d+\sqrt{a} e\right ) \log \left (\sqrt{a}-\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{c} x^2\right )}{16 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}-\frac{\left (3 \sqrt{c} d+\sqrt{a} e\right ) \log \left (\sqrt{a}+\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{c} x^2\right )}{16 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}+\frac{\left (d^2 \left (\frac{\sqrt{c} d}{\sqrt{a}}-e\right )\right ) \int \frac{1}{\frac{\sqrt{a}}{\sqrt{c}}-\frac{\sqrt{2} \sqrt [4]{a} x}{\sqrt [4]{c}}+x^2} \, dx}{4 \left (c d^2+a e^2\right )^2}+\frac{\left (d^2 \left (\frac{\sqrt{c} d}{\sqrt{a}}-e\right )\right ) \int \frac{1}{\frac{\sqrt{a}}{\sqrt{c}}+\frac{\sqrt{2} \sqrt [4]{a} x}{\sqrt [4]{c}}+x^2} \, dx}{4 \left (c d^2+a e^2\right )^2}-\frac{\left (\sqrt [4]{c} d^2 \left (\sqrt{c} d+\sqrt{a} e\right )\right ) \int \frac{\frac{\sqrt{2} \sqrt [4]{a}}{\sqrt [4]{c}}+2 x}{-\frac{\sqrt{a}}{\sqrt{c}}-\frac{\sqrt{2} \sqrt [4]{a} x}{\sqrt [4]{c}}-x^2} \, dx}{4 \sqrt{2} a^{3/4} \left (c d^2+a e^2\right )^2}-\frac{\left (\sqrt [4]{c} d^2 \left (\sqrt{c} d+\sqrt{a} e\right )\right ) \int \frac{\frac{\sqrt{2} \sqrt [4]{a}}{\sqrt [4]{c}}-2 x}{-\frac{\sqrt{a}}{\sqrt{c}}+\frac{\sqrt{2} \sqrt [4]{a} x}{\sqrt [4]{c}}-x^2} \, dx}{4 \sqrt{2} a^{3/4} \left (c d^2+a e^2\right )^2}-\frac{\left (3 \sqrt{c} d-\sqrt{a} e\right ) \operatorname{Subst}\left (\int \frac{1}{-1-x^2} \, dx,x,1-\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}+\frac{\left (3 \sqrt{c} d-\sqrt{a} e\right ) \operatorname{Subst}\left (\int \frac{1}{-1-x^2} \, dx,x,1+\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}\\ &=-\frac{x \left (d-e x^2\right )}{4 \left (c d^2+a e^2\right ) \left (a+c x^4\right )}+\frac{d^{3/2} e^{3/2} \tan ^{-1}\left (\frac{\sqrt{e} x}{\sqrt{d}}\right )}{\left (c d^2+a e^2\right )^2}+\frac{\left (3 \sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (1-\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}-\frac{\left (3 \sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (1+\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}-\frac{\sqrt [4]{c} d^2 \left (\sqrt{c} d+\sqrt{a} e\right ) \log \left (\sqrt{a}-\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{c} x^2\right )}{4 \sqrt{2} a^{3/4} \left (c d^2+a e^2\right )^2}+\frac{\left (3 \sqrt{c} d+\sqrt{a} e\right ) \log \left (\sqrt{a}-\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{c} x^2\right )}{16 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}+\frac{\sqrt [4]{c} d^2 \left (\sqrt{c} d+\sqrt{a} e\right ) \log \left (\sqrt{a}+\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{c} x^2\right )}{4 \sqrt{2} a^{3/4} \left (c d^2+a e^2\right )^2}-\frac{\left (3 \sqrt{c} d+\sqrt{a} e\right ) \log \left (\sqrt{a}+\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{c} x^2\right )}{16 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}+\frac{\left (\sqrt [4]{c} d^2 \left (\sqrt{c} d-\sqrt{a} e\right )\right ) \operatorname{Subst}\left (\int \frac{1}{-1-x^2} \, dx,x,1-\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt{2} a^{3/4} \left (c d^2+a e^2\right )^2}-\frac{\left (\sqrt [4]{c} d^2 \left (\sqrt{c} d-\sqrt{a} e\right )\right ) \operatorname{Subst}\left (\int \frac{1}{-1-x^2} \, dx,x,1+\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt{2} a^{3/4} \left (c d^2+a e^2\right )^2}\\ &=-\frac{x \left (d-e x^2\right )}{4 \left (c d^2+a e^2\right ) \left (a+c x^4\right )}+\frac{d^{3/2} e^{3/2} \tan ^{-1}\left (\frac{\sqrt{e} x}{\sqrt{d}}\right )}{\left (c d^2+a e^2\right )^2}-\frac{\sqrt [4]{c} d^2 \left (\sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (1-\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt{2} a^{3/4} \left (c d^2+a e^2\right )^2}+\frac{\left (3 \sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (1-\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}+\frac{\sqrt [4]{c} d^2 \left (\sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (1+\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt{2} a^{3/4} \left (c d^2+a e^2\right )^2}-\frac{\left (3 \sqrt{c} d-\sqrt{a} e\right ) \tan ^{-1}\left (1+\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}-\frac{\sqrt [4]{c} d^2 \left (\sqrt{c} d+\sqrt{a} e\right ) \log \left (\sqrt{a}-\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{c} x^2\right )}{4 \sqrt{2} a^{3/4} \left (c d^2+a e^2\right )^2}+\frac{\left (3 \sqrt{c} d+\sqrt{a} e\right ) \log \left (\sqrt{a}-\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{c} x^2\right )}{16 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}+\frac{\sqrt [4]{c} d^2 \left (\sqrt{c} d+\sqrt{a} e\right ) \log \left (\sqrt{a}+\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{c} x^2\right )}{4 \sqrt{2} a^{3/4} \left (c d^2+a e^2\right )^2}-\frac{\left (3 \sqrt{c} d+\sqrt{a} e\right ) \log \left (\sqrt{a}+\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{c} x^2\right )}{16 \sqrt{2} a^{3/4} c^{3/4} \left (c d^2+a e^2\right )}\\ \end{align*}

Mathematica [A]  time = 0.285044, size = 423, normalized size = 0.62 \[ \frac{\frac{\sqrt{2} \left (a^{3/2} e^3-3 \sqrt{a} c d^2 e+3 a \sqrt{c} d e^2-c^{3/2} d^3\right ) \log \left (-\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{a}+\sqrt{c} x^2\right )}{a^{3/4} c^{3/4}}+\frac{\sqrt{2} \left (-a^{3/2} e^3+3 \sqrt{a} c d^2 e-3 a \sqrt{c} d e^2+c^{3/2} d^3\right ) \log \left (\sqrt{2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt{a}+\sqrt{c} x^2\right )}{a^{3/4} c^{3/4}}-\frac{2 \sqrt{2} \left (a^{3/2} e^3-3 \sqrt{a} c d^2 e-3 a \sqrt{c} d e^2+c^{3/2} d^3\right ) \tan ^{-1}\left (1-\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{a^{3/4} c^{3/4}}+\frac{2 \sqrt{2} \left (a^{3/2} e^3-3 \sqrt{a} c d^2 e-3 a \sqrt{c} d e^2+c^{3/2} d^3\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt [4]{c} x}{\sqrt [4]{a}}+1\right )}{a^{3/4} c^{3/4}}+\frac{8 \left (e x^3-d x\right ) \left (a e^2+c d^2\right )}{a+c x^4}+32 d^{3/2} e^{3/2} \tan ^{-1}\left (\frac{\sqrt{e} x}{\sqrt{d}}\right )}{32 \left (a e^2+c d^2\right )^2} \]

Antiderivative was successfully verified.

[In]

Integrate[x^4/((d + e*x^2)*(a + c*x^4)^2),x]

[Out]

((8*(c*d^2 + a*e^2)*(-(d*x) + e*x^3))/(a + c*x^4) + 32*d^(3/2)*e^(3/2)*ArcTan[(Sqrt[e]*x)/Sqrt[d]] - (2*Sqrt[2
]*(c^(3/2)*d^3 - 3*Sqrt[a]*c*d^2*e - 3*a*Sqrt[c]*d*e^2 + a^(3/2)*e^3)*ArcTan[1 - (Sqrt[2]*c^(1/4)*x)/a^(1/4)])
/(a^(3/4)*c^(3/4)) + (2*Sqrt[2]*(c^(3/2)*d^3 - 3*Sqrt[a]*c*d^2*e - 3*a*Sqrt[c]*d*e^2 + a^(3/2)*e^3)*ArcTan[1 +
 (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/(a^(3/4)*c^(3/4)) + (Sqrt[2]*(-(c^(3/2)*d^3) - 3*Sqrt[a]*c*d^2*e + 3*a*Sqrt[c]*
d*e^2 + a^(3/2)*e^3)*Log[Sqrt[a] - Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2])/(a^(3/4)*c^(3/4)) + (Sqrt[2]*(c^(
3/2)*d^3 + 3*Sqrt[a]*c*d^2*e - 3*a*Sqrt[c]*d*e^2 - a^(3/2)*e^3)*Log[Sqrt[a] + Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt
[c]*x^2])/(a^(3/4)*c^(3/4)))/(32*(c*d^2 + a*e^2)^2)

________________________________________________________________________________________

Maple [A]  time = 0.016, size = 848, normalized size = 1.2 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(x^4/(e*x^2+d)/(c*x^4+a)^2,x)

[Out]

1/4/(a*e^2+c*d^2)^2/(c*x^4+a)*x^3*e^3*a+1/4/(a*e^2+c*d^2)^2/(c*x^4+a)*x^3*e*d^2*c-1/4/(a*e^2+c*d^2)^2/(c*x^4+a
)*x*e^2*d*a-1/4/(a*e^2+c*d^2)^2/(c*x^4+a)*x*c*d^3-3/16/(a*e^2+c*d^2)^2*(1/c*a)^(1/4)*2^(1/2)*arctan(2^(1/2)/(1
/c*a)^(1/4)*x-1)*d*e^2+1/16/(a*e^2+c*d^2)^2*(1/c*a)^(1/4)/a*2^(1/2)*arctan(2^(1/2)/(1/c*a)^(1/4)*x-1)*c*d^3-3/
32/(a*e^2+c*d^2)^2*(1/c*a)^(1/4)*2^(1/2)*ln((x^2+(1/c*a)^(1/4)*x*2^(1/2)+(1/c*a)^(1/2))/(x^2-(1/c*a)^(1/4)*x*2
^(1/2)+(1/c*a)^(1/2)))*d*e^2+1/32/(a*e^2+c*d^2)^2*(1/c*a)^(1/4)/a*2^(1/2)*ln((x^2+(1/c*a)^(1/4)*x*2^(1/2)+(1/c
*a)^(1/2))/(x^2-(1/c*a)^(1/4)*x*2^(1/2)+(1/c*a)^(1/2)))*c*d^3-3/16/(a*e^2+c*d^2)^2*(1/c*a)^(1/4)*2^(1/2)*arcta
n(2^(1/2)/(1/c*a)^(1/4)*x+1)*d*e^2+1/16/(a*e^2+c*d^2)^2*(1/c*a)^(1/4)/a*2^(1/2)*arctan(2^(1/2)/(1/c*a)^(1/4)*x
+1)*c*d^3+1/16/(a*e^2+c*d^2)^2/c/(1/c*a)^(1/4)*2^(1/2)*arctan(2^(1/2)/(1/c*a)^(1/4)*x-1)*a*e^3-3/16/(a*e^2+c*d
^2)^2/(1/c*a)^(1/4)*2^(1/2)*arctan(2^(1/2)/(1/c*a)^(1/4)*x-1)*d^2*e+1/32/(a*e^2+c*d^2)^2/c/(1/c*a)^(1/4)*2^(1/
2)*ln((x^2-(1/c*a)^(1/4)*x*2^(1/2)+(1/c*a)^(1/2))/(x^2+(1/c*a)^(1/4)*x*2^(1/2)+(1/c*a)^(1/2)))*a*e^3-3/32/(a*e
^2+c*d^2)^2/(1/c*a)^(1/4)*2^(1/2)*ln((x^2-(1/c*a)^(1/4)*x*2^(1/2)+(1/c*a)^(1/2))/(x^2+(1/c*a)^(1/4)*x*2^(1/2)+
(1/c*a)^(1/2)))*d^2*e+1/16/(a*e^2+c*d^2)^2/c/(1/c*a)^(1/4)*2^(1/2)*arctan(2^(1/2)/(1/c*a)^(1/4)*x+1)*a*e^3-3/1
6/(a*e^2+c*d^2)^2/(1/c*a)^(1/4)*2^(1/2)*arctan(2^(1/2)/(1/c*a)^(1/4)*x+1)*d^2*e+d^2*e^2/(a*e^2+c*d^2)^2/(d*e)^
(1/2)*arctan(e*x/(d*e)^(1/2))

________________________________________________________________________________________

Maxima [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: ValueError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^4/(e*x^2+d)/(c*x^4+a)^2,x, algorithm="maxima")

[Out]

Exception raised: ValueError

________________________________________________________________________________________

Fricas [B]  time = 59.1055, size = 19637, normalized size = 28.67 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^4/(e*x^2+d)/(c*x^4+a)^2,x, algorithm="fricas")

[Out]

[1/16*(4*(c*d^2*e + a*e^3)*x^3 - (a*c^2*d^4 + 2*a^2*c*d^2*e^2 + a^3*e^4 + (c^3*d^4 + 2*a*c^2*d^2*e^2 + a^2*c*e
^4)*x^4)*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 + (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4
 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*
e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^
12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2
*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8))*
log(-(c^4*d^8 - 14*a*c^3*d^6*e^2 + 14*a^3*c*d^2*e^6 - a^4*e^8)*x + (a*c^5*d^9 - 18*a^2*c^4*d^7*e^2 + 60*a^3*c^
3*d^5*e^4 - 46*a^4*c^2*d^3*e^6 + 3*a^5*c*d*e^8 + (3*a^3*c^7*d^10*e + 11*a^4*c^6*d^8*e^3 + 14*a^5*c^5*d^6*e^5 +
 6*a^6*c^4*d^4*e^7 - a^7*c^3*d^2*e^9 - a^8*c^2*e^11)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4
 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*
e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^
12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 + (a*c^5*d^8 + 4*
a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255
*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 +
8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 2
8*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4
 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8))) + (a*c^2*d^4 + 2*a^2*c*d^2*e^2 + a^3*e^4 + (c^3*d^4 + 2*a*c^2*d^2*e^2 + a^
2*c*e^4)*x^4)*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 + (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^
4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3
*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c
^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^
4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e
^8))*log(-(c^4*d^8 - 14*a*c^3*d^6*e^2 + 14*a^3*c*d^2*e^6 - a^4*e^8)*x - (a*c^5*d^9 - 18*a^2*c^4*d^7*e^2 + 60*a
^3*c^3*d^5*e^4 - 46*a^4*c^2*d^3*e^6 + 3*a^5*c*d*e^8 + (3*a^3*c^7*d^10*e + 11*a^4*c^6*d^8*e^3 + 14*a^5*c^5*d^6*
e^5 + 6*a^6*c^4*d^4*e^7 - a^7*c^3*d^2*e^9 - a^8*c^2*e^11)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^
8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*
d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d
^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 + (a*c^5*d^8
 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2
+ 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^
16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^1
0 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^
4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8))) - (a*c^2*d^4 + 2*a^2*c*d^2*e^2 + a^3*e^4 + (c^3*d^4 + 2*a*c^2*d^2*e^2
 + a^2*c*e^4)*x^4)*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 - (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c
^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^
3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*
a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^
10*c^4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^
5*c*e^8))*log(-(c^4*d^8 - 14*a*c^3*d^6*e^2 + 14*a^3*c*d^2*e^6 - a^4*e^8)*x + (a*c^5*d^9 - 18*a^2*c^4*d^7*e^2 +
 60*a^3*c^3*d^5*e^4 - 46*a^4*c^2*d^3*e^6 + 3*a^5*c*d*e^8 - (3*a^3*c^7*d^10*e + 11*a^4*c^6*d^8*e^3 + 14*a^5*c^5
*d^6*e^5 + 6*a^6*c^4*d^4*e^7 - a^7*c^3*d^2*e^9 - a^8*c^2*e^11)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c
^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*
c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*
c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 - (a*c^
5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10
*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^
11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^
6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c
^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8))) + (a*c^2*d^4 + 2*a^2*c*d^2*e^2 + a^3*e^4 + (c^3*d^4 + 2*a*c^2*d^
2*e^2 + a^2*c*e^4)*x^4)*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 - (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*
a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 4
52*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2
+ 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 +
 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6
 + a^5*c*e^8))*log(-(c^4*d^8 - 14*a*c^3*d^6*e^2 + 14*a^3*c*d^2*e^6 - a^4*e^8)*x - (a*c^5*d^9 - 18*a^2*c^4*d^7*
e^2 + 60*a^3*c^3*d^5*e^4 - 46*a^4*c^2*d^3*e^6 + 3*a^5*c*d*e^8 - (3*a^3*c^7*d^10*e + 11*a^4*c^6*d^8*e^3 + 14*a^
5*c^5*d^6*e^5 + 6*a^6*c^4*d^4*e^7 - a^7*c^3*d^2*e^9 - a^8*c^2*e^11)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*
a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8
*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28
*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 -
(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5
*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a
^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c
^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*
a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8))) + 8*(c*d*e*x^4 + a*d*e)*sqrt(-d*e)*log((e*x^2 + 2*sqrt(-d*e
)*x - d)/(e*x^2 + d)) - 4*(c*d^3 + a*d*e^2)*x)/(a*c^2*d^4 + 2*a^2*c*d^2*e^2 + a^3*e^4 + (c^3*d^4 + 2*a*c^2*d^2
*e^2 + a^2*c*e^4)*x^4), 1/16*(4*(c*d^2*e + a*e^3)*x^3 + 16*(c*d*e*x^4 + a*d*e)*sqrt(d*e)*arctan(sqrt(d*e)*x/d)
 - (a*c^2*d^4 + 2*a^2*c*d^2*e^2 + a^3*e^4 + (c^3*d^4 + 2*a*c^2*d^2*e^2 + a^2*c*e^4)*x^4)*sqrt((6*c^2*d^5*e - 2
0*a*c*d^3*e^3 + 6*a^2*d*e^5 + (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e
^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30
*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 +
 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^
5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8))*log(-(c^4*d^8 - 14*a*c^3*d^6*e
^2 + 14*a^3*c*d^2*e^6 - a^4*e^8)*x + (a*c^5*d^9 - 18*a^2*c^4*d^7*e^2 + 60*a^3*c^3*d^5*e^4 - 46*a^4*c^2*d^3*e^6
 + 3*a^5*c*d*e^8 + (3*a^3*c^7*d^10*e + 11*a^4*c^6*d^8*e^3 + 14*a^5*c^5*d^6*e^5 + 6*a^6*c^4*d^4*e^7 - a^7*c^3*d
^2*e^9 - a^8*c^2*e^11)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a
^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56
*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^1
1*c^3*e^16)))*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 + (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^
4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3
*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c
^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^
4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e
^8))) + (a*c^2*d^4 + 2*a^2*c*d^2*e^2 + a^3*e^4 + (c^3*d^4 + 2*a*c^2*d^2*e^2 + a^2*c*e^4)*x^4)*sqrt((6*c^2*d^5*
e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 + (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^
5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8
 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*
e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))/
(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8))*log(-(c^4*d^8 - 14*a*c^3*
d^6*e^2 + 14*a^3*c*d^2*e^6 - a^4*e^8)*x - (a*c^5*d^9 - 18*a^2*c^4*d^7*e^2 + 60*a^3*c^3*d^5*e^4 - 46*a^4*c^2*d^
3*e^6 + 3*a^5*c*d*e^8 + (3*a^3*c^7*d^10*e + 11*a^4*c^6*d^8*e^3 + 14*a^5*c^5*d^6*e^5 + 6*a^6*c^4*d^4*e^7 - a^7*
c^3*d^2*e^9 - a^8*c^2*e^11)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 +
255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4
 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14
+ a^11*c^3*e^16)))*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 + (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c
^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^
3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*
a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^
10*c^4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^
5*c*e^8))) - (a*c^2*d^4 + 2*a^2*c*d^2*e^2 + a^3*e^4 + (c^3*d^4 + 2*a*c^2*d^2*e^2 + a^2*c*e^4)*x^4)*sqrt((6*c^2
*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 - (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6
 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^
4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*
d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^1
6)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8))*log(-(c^4*d^8 - 14*a
*c^3*d^6*e^2 + 14*a^3*c*d^2*e^6 - a^4*e^8)*x + (a*c^5*d^9 - 18*a^2*c^4*d^7*e^2 + 60*a^3*c^3*d^5*e^4 - 46*a^4*c
^2*d^3*e^6 + 3*a^5*c*d*e^8 - (3*a^3*c^7*d^10*e + 11*a^4*c^6*d^8*e^3 + 14*a^5*c^5*d^6*e^5 + 6*a^6*c^4*d^4*e^7 -
 a^7*c^3*d^2*e^9 - a^8*c^2*e^11)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e
^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^1
2*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*
e^14 + a^11*c^3*e^16)))*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 - (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*
a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 4
52*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2
+ 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 +
 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6
 + a^5*c*e^8))) + (a*c^2*d^4 + 2*a^2*c*d^2*e^2 + a^3*e^4 + (c^3*d^4 + 2*a*c^2*d^2*e^2 + a^2*c*e^4)*x^4)*sqrt((
6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 - (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^
2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c
^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6
*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^
3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8))*log(-(c^4*d^8 -
 14*a*c^3*d^6*e^2 + 14*a^3*c*d^2*e^6 - a^4*e^8)*x - (a*c^5*d^9 - 18*a^2*c^4*d^7*e^2 + 60*a^3*c^3*d^5*e^4 - 46*
a^4*c^2*d^3*e^6 + 3*a^5*c*d*e^8 - (3*a^3*c^7*d^10*e + 11*a^4*c^6*d^8*e^3 + 14*a^5*c^5*d^6*e^5 + 6*a^6*c^4*d^4*
e^7 - a^7*c^3*d^2*e^9 - a^8*c^2*e^11)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^4 - 452*a^3*c^3*
d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14*e^2 + 28*a^5*c^
9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e^12 + 8*a^10*c^4
*d^2*e^14 + a^11*c^3*e^16)))*sqrt((6*c^2*d^5*e - 20*a*c*d^3*e^3 + 6*a^2*d*e^5 - (a*c^5*d^8 + 4*a^2*c^4*d^6*e^2
 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)*sqrt(-(c^6*d^12 - 30*a*c^5*d^10*e^2 + 255*a^2*c^4*d^8*e^
4 - 452*a^3*c^3*d^6*e^6 + 255*a^4*c^2*d^4*e^8 - 30*a^5*c*d^2*e^10 + a^6*e^12)/(a^3*c^11*d^16 + 8*a^4*c^10*d^14
*e^2 + 28*a^5*c^9*d^12*e^4 + 56*a^6*c^8*d^10*e^6 + 70*a^7*c^7*d^8*e^8 + 56*a^8*c^6*d^6*e^10 + 28*a^9*c^5*d^4*e
^12 + 8*a^10*c^4*d^2*e^14 + a^11*c^3*e^16)))/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^
2*e^6 + a^5*c*e^8))) - 4*(c*d^3 + a*d*e^2)*x)/(a*c^2*d^4 + 2*a^2*c*d^2*e^2 + a^3*e^4 + (c^3*d^4 + 2*a*c^2*d^2*
e^2 + a^2*c*e^4)*x^4)]

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x**4/(e*x**2+d)/(c*x**4+a)**2,x)

[Out]

Timed out

________________________________________________________________________________________

Giac [A]  time = 1.15073, size = 791, normalized size = 1.15 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^4/(e*x^2+d)/(c*x^4+a)^2,x, algorithm="giac")

[Out]

d^(3/2)*arctan(x*e^(1/2)/sqrt(d))*e^(3/2)/(c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4) + 1/8*((a*c^3)^(1/4)*c^3*d^3 - 3
*(a*c^3)^(1/4)*a*c^2*d*e^2 - 3*(a*c^3)^(3/4)*c*d^2*e + (a*c^3)^(3/4)*a*e^3)*arctan(1/2*sqrt(2)*(2*x + sqrt(2)*
(a/c)^(1/4))/(a/c)^(1/4))/(sqrt(2)*a*c^5*d^4 + 2*sqrt(2)*a^2*c^4*d^2*e^2 + sqrt(2)*a^3*c^3*e^4) + 1/8*((a*c^3)
^(1/4)*c^3*d^3 - 3*(a*c^3)^(1/4)*a*c^2*d*e^2 - 3*(a*c^3)^(3/4)*c*d^2*e + (a*c^3)^(3/4)*a*e^3)*arctan(1/2*sqrt(
2)*(2*x - sqrt(2)*(a/c)^(1/4))/(a/c)^(1/4))/(sqrt(2)*a*c^5*d^4 + 2*sqrt(2)*a^2*c^4*d^2*e^2 + sqrt(2)*a^3*c^3*e
^4) + 1/16*((a*c^3)^(1/4)*c^3*d^3 - 3*(a*c^3)^(1/4)*a*c^2*d*e^2 + 3*(a*c^3)^(3/4)*c*d^2*e - (a*c^3)^(3/4)*a*e^
3)*log(x^2 + sqrt(2)*x*(a/c)^(1/4) + sqrt(a/c))/(sqrt(2)*a*c^5*d^4 + 2*sqrt(2)*a^2*c^4*d^2*e^2 + sqrt(2)*a^3*c
^3*e^4) - 1/16*((a*c^3)^(1/4)*c^3*d^3 - 3*(a*c^3)^(1/4)*a*c^2*d*e^2 + 3*(a*c^3)^(3/4)*c*d^2*e - (a*c^3)^(3/4)*
a*e^3)*log(x^2 - sqrt(2)*x*(a/c)^(1/4) + sqrt(a/c))/(sqrt(2)*a*c^5*d^4 + 2*sqrt(2)*a^2*c^4*d^2*e^2 + sqrt(2)*a
^3*c^3*e^4) + 1/4*(x^3*e - d*x)/((c*x^4 + a)*(c*d^2 + a*e^2))